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The sensitivity of inertial confinement fusion implosions of the type performed on the National 
Ignition Facility (NIF) HI to low-mode flux asymmetries has been investigated numerically. It 
is shown that large-amplitude, low-order mode shapes (Legendre polynomial P4), resulting from 
associated low order flux asymmetries, cause spatial variations in capsule & fuel momentum that 
prevent the DT "ice" layer from being decelerated uniformly by the hot spot pressure. This reduces 
the transfer of kinetic to internal energy of the central hot spot, thus reducing neutron yield. 
Furthermore, synthetic gated x-ray images of the hot spot self-emission indicate that P4 shapes may 
be unquantifiable for DT layered capsules. Instead the positive P4 asymmetry "aliases" itself as an 
oblate P2 in the x-ray self emission images. Correction of this apparent P2 distortion can further 
distort the implosion while creating a round x-ray image. Long wavelength asymmetries may be 
playing a significant role in the observed yield reduction of NIF DT implosions relative to detailed 
post-shot 2D simulations. 



Indirect-drive inertial confinement fusion (ICF) [T]-[3] 
uses lasers to heat the inside of a cavity (or hohlraum). 
The absorbed laser energy is re-emitted as approximately 
black-body radiation in the soft x-ray regime. These x- 
rays heat the outer surface of a hollow, spherical, ~ 2 
mm diameter, shell that contains a ~ 70 pm thick layer 
of cryogenically frozen Deuterium and Tritium fuel ( "DT 
fuel" or "DT layered capsules"). The heated outer shell 
ablates, which creates a reaction force, accelerating the 
remaining shell spherically inwards at extremely high ve- 
locity (^ 350 km/s). During the implosion, spherical 
convergence causes the pressure in the gaseous void (or 
hot spot) within the shell to rise rapidly. This pressure 
decelerates the shell, simultaneously compressing the 
solid fuel and converting the shell's kinetic energy into 
hot spot internal energy. If this conversion rate exceeds 
loss rates due to thermal conduction and bremsstrahlung 
radiation, the hot spot will heat, initiating DT fusion re- 
actions. Provided the hot spot areal density is sufficient, 
a-particles created by the fusion reactions will redeposit 
their energy locally, further heating the hot spot, result- 
ing in bootstrap heating, ignition, and propagation of 
burn into the surrounding cold fuel. Numerical modeling 
indicates that the National Ignition Facility (NIF) can, 
for the first time, initiate inertial fusion ignition in the 
laboratory [41^. 

In this Letter, the effects of large, low-mode asymme- 
tries in the x-ray drive are examined numerically. The 
non-uniformity of the x-ray flux incident upon the shell 
and the resultant non-spherical shell shapes can be de- 
scribed mathematically as a series of Legendre polynomi- 
als [7]. It is shown that a large- amplitude P4 implosion 
asymmetry, that might result from low-order hohlraum 



generated fiux asymmetries, causes spatial variations in 
the capsule & fuel momentum. This can inhibit the DT 
fuel from being decelerated uniformly by the hot spot 
pressure, reducing the efficient transfer of implosion ki- 
netic energy to hot spot internal energy thus significantly 
reducing the capsule performance. Furthermore, simu- 
lated gated x-ray images of the hot spot self-emission 
show reduced sensitivity to the P4 mode, instead the im- 
ages appear to have a pronounced oblate P2 shape. Re- 
ducing the amplitude of the oblate P2 shape (as measured 
from the x-ray image) further reduces the sensitivity to 
the P4 mode such that no quantitative evaluation of the 
hot spot a4 (where a4 is the amplitude of the P4 mode) 
can be made, furthermore the x-ray images are circular 
despite the capsule shape being highly distorted. Com- 
parisons are made between key physical properties of the 
implosion, synthetically generated experimental observ- 
ables, and NIF experimental data. 

The indirect-drive approach to ICF smooths high mode 
spatial non-uniformities in the x-ray flux incident on the 
capsule, however the spatial distribution of the cones 
of laser beams which illuminate the hohlraum means 
that low mode x-ray flux non-uniformities can occur [T], 
these are considerably lower mode than those recently 
examined by Thomas et al [8]. The growth of Legen- 
dre polynomial P4 capsule shapes was investigated using 
the radiation-magnetohydrodynamics code Hydra [9 . A 
frequency and time dependent x-ray source was devel- 
oped to drive these capsule-only simulations. The ini- 
tial x-ray drive was taken from an integrated hohlraum 
simulation and then adjusted to match the shock timing 
data obtained using the VISAR diagnostic [10] from NIF 
shot Nl 10521, and the capsule implosion trajectory [11] 
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Figure 1. (a) Applied radiation temperature vs time. Colors 
depict when during the drive pulse the various flux asymme- 
tries were applied as shown on subsequent plots, (b) The 
relationship between applied flux asymmetry and hot spot a4 
at x-ray bangtime (the time at which peak emission occurs). 
Each point represents a 2D Hydra run. (c) Burn averaged hot 
spot + fuel + ablator pr vs hot spot a4 at x-ray bangtime: 
large dots are spatially averaged pr, while the smaller points 
with the same color and a4 are the maxa and minima of the 
spatially averaged value; large variations in pr occur due to 
P4. (d) The burn averaged energy partition as a function of 
hot spot a4; increasing P4 perturbations prevent kinetic en- 
ergy (black) from being converted to both hot spot (red) and 
solid fuel (blue) internal energy during stagnation, (e) Burn 
averaged hot spot pressure as a function of hot spot a^. (f) 
Total thermonuclear neutron yield as a function of hot spot 
a4; yield varies by a factor of 15. 



measured on NIF shot N110625. A 90° 'wedge' of the 
capsule (27r Sr) was modeled using two-dimensional (2D) 
cylindrically-symmetric geometry with 256 x 312 cells. 
Doubling and quadrupling the cell resolution demon- 
strated convergence. In all runs the Quotidian Equation 



of State [12] was used with tabular opacities and multi- 
group radiation diffusion. The effects of Legendre poly- 
nomial P4 hohlraum flux asymmetries were investigated 
by perturbing the tuned x-ray drive with spatially vary- 
ing flux asymmetries of the form: fdsiO^t) = (aoPo + 
a4P4{0)) * fdsit) where fds is the energy density of the 
tuned photon frequency dependent x-ray drive source, 
is the amplitude of the n^^ Legendre polynomial, ao = 1, 
a4 = (±0.10, ±0.05, ±0.025, ±0.01, ±0.005, ±0.0025), 
the angle between the equatorial plane and polar axis, 
and t time. Hydra modeling of the hohlraum & capsule 
for nominal implosions suggests the flux asymmetry inci- 
dent on the capsule would be expected to vary by < 3% 
except for in the first ~ 2 ns of the laser pulse were it can 
vary by up to 10% [13 . The flux asymmetries were ap- 
plied 100 pm from the capsule ablation front during dis- 
crete time intervals (see fig. [ija)), creating a database 
of > 200 2D modeling runs of both DT layered implo- 
sions and DHe^ gas filled capsules with a surrogate inner 
CH layer of equal mass to a DT fuel layer (symmetry cap- 
sules). Time resolved synthetic gated x-ray images of the 
hot spot self-emission > 6 keV, including its attenuation 
by the compressed fuel and ablator, were created from 
both polar and equatorial directions by post processing 
each Hydra run. The images were blurred in order to 
reproduce the 11 pm resolution of the diagnostic. The 
key implosion performance metrics (neutron yield, hot 
spot pressure, mass, volume, density, ion & electron tem- 
peratures, the effective ion temperature computed from 
the FWHM of the DT neutron spectrum, fuel and abla- 
tor areal density (pr) and kinetic energy) were extracted 
from the simulations. The hot spot shape was evaluated 
as a function of time by performing a Legendre poly- 
nomial decomposition (modes 1-10) of the appropriate 
contour. For DT layered capsules the hot spot contour is 
defined for each angular 'strip' of cells j as the minimum 
radius where T^. > ^T^.^^^ and pj < \pj^^^ where Tg is 
the electron temperature and p the mass density, 'max' 
denotes the maximum value within the j^^ strip. This 
has been found to produce a robust definition of the hot 
spot even for highly distorted implosions. The 17% con- 
tour of the gated x-ray diagnostic (GXD) is used both for 
the synthetic GXD and experimentally, as previous stud- 
ies have shown this provides a faithful representation of 
the hot spot shape for small departures from sphericity. 

The applied Legendre P4 flux asymmetries induce P4 
hot spot shapes at stagnation (see Figs. [2] (a) & (c)), 
the sign of which is dependent on the timing of the ap- 
plied flux asymmetry. If the asymmetry is present only 
during the shock compression phase (the first ~ 18 ns), 
shocks created in regions of the capsule exposed to higher 
flux propagate faster, these faster shocks break out of the 
inner DT ice layer earlier, causing these regions to move 
ahead of those exposed to less flux. This also causes abla- 
tor mass to flow laterally, away from the high flux region. 
Consequently during peak drive the regions initially ex- 
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Figure 3. (a) DT layered capsule hot spot 04 plotted against 
the synthetic GXD 04; particularly for large positive a4 the 
GXD is unable to effectively measure the amplitude of the P4 
mode, (b) DT layered capsule hot spot a4 plotted against the 
synthetic GXD a2; the GXD measures a significant P2 mode 
amplitude despite the DT layered capsule hot spot a2 being 
lb 1 Jim (not shown) . 
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Figure 2. Axis of rotational symmetry is horizontal at Radius 
= pm. (a) DT layered capsule density plot at x-ray bang- 
time showing a positive Legendre polynomial P4 shape. This 
simulation had a 10% flux asymmetry applied from 11.5-14 
ns. Black arrows indicate the mass flows which occur dur- 
ing stagnation. After bangtime 'fingers' of fuel continue to 
flow inwards (red arrows). White dots depict the hot spot 
contour, (b) Synthetic gated x-ray image of the hot spot self 
emission from (a), white dots show the 17% contour, (c) Fig. 
(a) 100 ps later: due to burn truncation with large this is 
the neutron bangtime for an equivalent spherical implosion, 
(d) The synthetic GXD from (c), showing a large negative P2 
and almost zero a4 despite the obvious P4 in (c). 

posed to high flux are at smaller radii, meaning they are 
accelerated less efficiently by the hohlraum flux and gain 
less total momentum. They can also have less ablator 
aerial density. The net effect is that the regions experi- 
encing high flux during shock compression will protrude 
outwards at stagnation. Conversely if the flux asymme- 
try is applied during peak drive, the regions of the capsule 
exposed to more flux gain more momentum, and pro- 
trude inwards at stagnation. Regardless of the timing of 
the applied asymmetry, during the stagnation phase of 
the implosion, pressure within the lower density hot spot 
decelerates the higher density fuel from peak velocity, 
making this interface Ray leigh- Taylor unstable [HI [15] . 
The instability will accentuate any shape imperfections 
during deceleration, as indicated by the significant simu- 
lated growth shown in figs. |2|a) & (c). 

The scalings of some important DT layered capsule 
implosion parameters as a function of hot spot mea- 
sured at x-ray bangtime are summarized in figure [l] 
Fig. [ijb) shows the relationship between applied P4 
flux perturbation amplitude and the resulting shape 
a4 at x-ray bangtime. Fig. [TJc) depicts the 'burn 



averaged' pr (the burn average of a quantity = 
(X]t=o^ QtEprodrdt) / j^Z^ Eprodr dt where Qt is Q at 
time t and Eprodr the thermonuclear energy production 
rate in time dt) as a function of hot spot a/^. Although 
the spatially averaged pr is relatively constant, the lat- 
eral mass flows caused by the P4 can create large spatial 
variations in pr. The regions with higher momentum 
continue to propagate radially inwards; fig. [TJd) depicts 
the remaining capsule kinetic energy as a function of a4 
and the partition of that energy into hot spot and solid 
fuel internal energy. For large a/^ less of the implosion 
kinetic energy is converted into hot spot internal energy 
and the hot spot pressure is reduced (see fig. [ije)). The 
reduction in neutron yield can be as large as 15x for hot 
spot a4 = 20 pm (flux asymmetry ~ 10%) as shown in 

fig- W))- 

Analysis of synthetic GXD images created from the 
2D Hydra runs suggest that the a4 measured experi- 
mentally with the GXD is not a true representation of 
the hot spot a4, particularly for large positive 04 ampli- 
tudes. Fig. [sja) depicts the relationship between the 
DT layered capsule "hot spot a4" (as previously defined) 
and that from the 17% contour of the synthetic GXD 
(the synthetic GXD a4), both were extracted at x-ray 
bangtime (the principal value used for analysis of exper- 
imental data). a4 measured from the synthetic GXD is 
consistently lower than that of the hot spot. The insensi- 
tivity to positive hot spot a4 is caused by lateral ablator 
mass flows which accumulate at ~ 45° at the expense 
of ablator material near the equator and poles (see Fig. 
[2] (a)). The ablator material is rotationally symmetric 
about the horizontal axis, so the accumulated material 
absorbs the x-rays emitted from the polar-lobes of the 
hot spot (left and right), while allowing x-rays to more 
readily pass through the equatorial regions (top and bot- 
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torn) . Consequently the polar-lobes of the hot spot which 
are visible in the density plots of Fig. |2] as dark regions 
(the hot spot is the central region of low density) are al- 
most completely invisible in the GXD plots compared to 
the emission through the equator. This causes the x-ray 
image to have a pronounced negative P2 shape (oblate 
or "pancaked"). As the hot spot a2 = =b 1 pm (a2 
is the amplitude of the P2 mode) for all these pure P4 
modelling runs, the P2 inferred from the x-ray image is a 
"false" negative P2 mode. This suggests that a negative 
P2 mode measured from the self-emission x-ray image 
may in fact be a signature of a positive P4 mode, al- 
though it does not, of course, preclude the presence of a 
true P2 mode. This is potentially important for interpre- 
tation of x-ray images from DT implosions, which often 
exhibit oblate (negative) P2 modes [16]. 

NIF experiments also use low convergence, DHe^ gas 
filled "symmetry capsules" which have a surrogate CH 
fuel mass. Using Hydra, symmetry capsule & DT lay- 
ered capsule pairs of runs were created by applying iden- 
tical x-ray drives to both capsules. When realistic exper- 
imental noise is applied to the synthetic GXD images, 
symmetry capsules, which have better signal to noise ra- 
tios and larger stagnation diameters, enable a far better 
measurement of the P4 mode than equivalent DT lay- 
ered capsules. Nevertheless, these calculations indicate 
symmetry capsules also show reduced sensitivity to a4, 
exhibit a "false" P2, and are quantitatively very similar 
to those shown in Figs, [s] (a) & (b) respectively. 

DT implosions on the NIF currently have yields 
3 — lOx below detailed 2D post-shot Hydra simulations 
[17] that match the measured shock timing, implosion 
velocity, and capsule and ice surface roughnesses. In 
comparison to experimentally measured or inferred val- 
ues [18] , these simulations have similar hot spot tempera- 
tures, however the hot spot volumes are reduced while the 
hot spot mass is increased, causing a 2 — 3x discrepancy 
in the hotspot density. P4 shape perturbations offer one 
mechanism which may explain these experimental obser- 
vations in particular bringing the yield and ion tempera- 
ture relationship into better agreement. In these simula- 
tions, the DT fuel and hot spot do not mix; clear bound- 
aries still exist (note these simulations use smooth cap- 
sules, but when nominal realistic capsule surface rough- 
ness [? ] was employed and modes up to 200 resolved, 
no significant implosion degradation occurred for the full 
range of a4). Consequently unlike high mode 'mix' [l] 
(where the hot spot can be radiatively cooled by high Z 
impurities), the simulated ion temperature inferred from 
the neutron spectrum remains unaffected at 3.9 ± 0.05 
keV for all a4. The large a/^ does however truncate the 
thermonuclear burn, moving both the neutron and x-ray 
bangtimes earlier in time, therefore as the capsule is still 
converging at bangtime, the hot spot size and volume are 
increased. The hot spot mass decreases with positive a4, 
bringing Hydra simulations approximately in line with 
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Figure 4. (a) Density plot of a DT layered capsule run with 
both P2 and P4 flux modes applied. Axis of rotational sym- 
metry is horizontal at Radius = pm. (b) The equatorial 
synthetic GXD image of (a) at the same time, dotted line 
shows 17% contour. Despite the highly non-spherical density 
distribution, the equatorial GXD image is almost perfectly 
round, (c) As the P2 flux amplitude is increased in order to 
make the GXD image look round (reducing 02) the sensitivity 
to a A is reduced towards zero. 



experimental data, as shown in Table |l] This compares 
NIF experimental data |18J with two Hydra implosions; 
one is perfectly spherical while the other has a hot spot 
a4 of +20 pm. Notable features are the significantly re- 
duced yield, reduced pressure, reduced hot spot mass, un- 
changed ion temperature and increased hot spot volume. 
We must emphasize, however, that this should not be in- 
terpreted as conclusive evidence that a P4 asymmetry is 
responsible for the observed reduced NIF capsule perfor- 
mance. Although this study has concentrated on the P4 
mode, it is likely that all low modes would reduce the 
conversion of capsule kinetic energy into hot spot pres- 
sure, and may result in similar ambiguity in the shape of 
the x-ray emission from the hot spot [19 . To explore the 
issue of low mode asymmetries further, experiments us- 
ing x-ray backlighters are currently being conducted on 
NIF to measure the implosion shape both in-flight pT] 
and at stagnation using Compton radiography [20j . 

As discussed, implosions with a significant P4 asym- 
metry can have a very apparent but "false" P2 asymme- 
try in GXD images. We find that attempting to correct 
this "false" P2 by increasing laser power to the hohlraum 
waist (capsule equator) can lead to a round GXD im- 
age even though the correction actually produces a more 
distorted DT fuel ice layer. This is depicted in fig. [4] 
for the case of a DT layered capsule where we applied 
and empirically adjusted a P2 flux asymmetry, in addi- 
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NIF expt. 


Hydra 


Hydra 




range 18 


(a4 = pm) 


(a4 = 20 pm) 


Pressure (GBar) 


57-81 


348 


115 


Mass (pg) 
Density (gcm~'^) 


2-6.4 


8 


5.5 


22-35 


136 


69 


Volume {xlO~^ cm^) 


0.9-1.9 


0.6 


1.0 


Tion (keV) 
Fuel pr (gcm~^) 


3.3-4.4 


3.9 


3.9 


0.77-0.98 


0.7 


0.72 


Yield (neut. xlO^^) 


1.9-6.0 


74 


5.3 



Table I. A comparison of NIF DT layered capsule experimen- 
tal data from 4 shots N110608-N1 10908 with two Hydra im- 
plosions, one spherical (a4 = pm, and another with a4 = 20 
pm. Large positive P4 brings the modeled implosion observ- 
ables approximately in line with the experimental data. 



the neutron yield by up to 15 x. Simulated x-ray images 
of the hot spot self-emission show reduced sensitivity to 
the positive P4 mode, instead the images appear to have 
a pronounced oblate P2 shape. Attempting to correct for 
this apparent P2 distortion can further distort the implo- 
sion while creating x-ray images which appear round & 
self-consistent from both equatorial and polar directions. 
This also further reduces the sensitivity to the P4 mode 
such that that no quantitative evaluation of the hot spot 
a4 can be made. Long wavelength asymmetries may be 
playing a significant role in the observed yield reduction 
of NIF DT implosions relative to detailed post-shot 2D 
simulations. 



tion to the original P4, in order to make the synthetic 
GXD image appear round. Fig. quantifies a related 
effect; as the applied P2 flux is increased in order to re- 
duce the "false" GXD a2 towards zero, there is a marked 
additional reduction in sensitivity to (relative to that 
shown in Fig. |3|. This suggests that attempts to tune 
the hohlraum to eliminate a "false" P2 can have the un- 
intended consequence of exacerbating overall asymmetry. 
Other information, such as comparison of the widths of 
images taken from both the polar and equatorial lines 
of sight [21 need to be taken into consideration. These 
simulations show that when a hotspot has a positive but 
pure P4 asymmetry the equatorial image width is larger 
than the width in the polar image (for negative P4 this 
is reversed). This could be used to identify an implosion 
where the measured P2 may be caused by a dominant P4 
asymmetry. However, our simulations also show us that 
the empirically but incorrectly tuned implosion of fig. 
[4] would have a polar image width that is equal to the 
equatorial image width, further misleading us into think- 
ing that we had engineered an approximately spherical 
implosion. A corollary of figure [4j is that it is possible to 
create imploded configurations which appear to be sym- 
metric in the GXD but, in fact, are significantly asym- 
metric and have greatly reduced performance in compar- 
ison to equivalent spherical implosions because a large 
fraction of the imploding shell's kinetic energy remains 
unstagnated. 

In summary, numerical simulations have been used to 
examine the sensitivity of implosions similar to those 
currently taking place on NIF to low-mode flux asym- 
metries. It is shown that Legendre polynomial P4 flux 
modes induce P4 shape modes at the time of capsule 
stagnation. The largest P4 amplitudes can cause up to 
50% of the capsule kinetic energy to remain unconverted 
to hot spot and DT ice internal energy, in turn reducing 
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